Influence of static electric fields on an optical ion trap 
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We recently reported on a proof-of-principle experiment demonstrating optical trapping of an ion 
in a single-beam dipole trap superimposed by a static electric potential [1]. Here, we first discuss the 
experimental procedures focussing on the influence and consequences of the static electric potential. 
These potentials can easily prevent successful optical trapping, if their configuration is not chosen 
carefully. Afterwards, we analyse the dipole trap experiments with different analytic models, in 
which different approximations are applied. According to these models the experimental results 
agree with recoil heating as the relevant heating effect. In addition, a Monte-Carlo simulation has 
been developed to refine the analysis. It reveals a large impact of the static electric potential on 
the dipole trap experiments in general. While it supports the results of the analytic models for the 
parameters used in the experiments, the analytic models cease their validity for significantly different 
parameters. Finally, we propose technical improvements for future realizations of experiments with 
optically trapped ions. 

PACS numbers: 37.10.-x 



INTRODUCTION 

Since the first trapping of ions in RF traps by Paul 
et al. in 1958 [2] and the first trapping of atoms in opti- 
cal traps by Chu et al. in 1986 [3], these trapping tech- 
niques have substantially advanced various fields of re- 
search. This includes, for example, quantum information 
processing [4], quantum metrology [5, 6], cold chemistry 
[7], and quantum degenerate gases [8]. However, trapping 
ions optically has been reported only recently [1]. A rea- 
son is that the Coulomb interaction leads to a high sen- 
sitivity to stray electric fields which might also account 
for the prejudice that ions would not be trappable in 
comparatively weak optical potentials. Potential depths 
of optical dipole traps, for example, are on the order of 
10 -3 K, while conventional RF traps have depths on the 
order of 10 4 K. 

An overview of the setup for optical trapping [1] is 
shown in Fig. 1. The basic experimental sequence is the 
following: Initially, a single ion is loaded and laser-cooled 
in a conventional, linear RF trap. Subsequently, a dipole 
trap beam is superimposed with the ion and the RF po- 
tential is ramped down to zero. Finally, after the dipole 
trap duration T optica i, the RF potential is ramped up 
again and the ion is detected in case of successful for- 
mer optical trapping. During all steps a static electric 
potential is retained, such that the total potential con- 
sists of the dipole trap and a static electric potential. 
Firstly, the static electric potential is required to com- 
pensate stray electric fields and, secondly, it prevents the 
ion from leaving the dipole trap in the propagation direc- 
tion of the beam k. Its focussing effect along the z-axis 
comes at the price of a defocussing effect in at least one 
radial direction, which has to be overcome by the dipole 
trap potential. Furthermore, the 45 ° angle between the 
z-axis and k results in an asymmetric total potential. 

In this paper we first discuss the experimental proce- 
dures and peculiarities. While the RF potential is still 
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Figure 1. Scheme of the trapping setup (color online). The 
linear RF trap consists of gold-plated, blade-like electrodes, 
which are mounted using alumina isolators and stainless steel 
holders. The DC electrodes are divided into six segments 
each: three belong to a zone with greater electrode-ion dis- 
tance and three to a zone with smaller electrode-ion distance. 
(The latter zone is not used in the current dipole trap exper- 
iments.) The atomic oven can evaporate a beam of neutral 
magnesium atoms mainly in direction of the zone with large 
electrode-ion distance. One or more atoms are photo-ionized 
and Doppler cooled in the RF potential using laser beams 
in the j/-z-plane, which intersect the z-axis at an angle of 
22.5 °. The dipole trap beam propagates in the a;-z-plane and 
crosses the z-axis at an angle of 45 ° . The fluorescence light 
from the ion(s) is detected with a CCD camera through a 
sapphire window above the trap (not shown). Parts of an iso- 
lator and a holder are omitted to offer a view on the trapping 
zone. The arrows of the coordinate system indicate a length 
scale of 10 mm; the distance between the atomic oven and its 
aperture is not to scale. 



on, the axial potential shows strong, unexpected anhar- 
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monicities. Using ions as sensors we measure, quantify 
and discuss the anharmonicities and deduce consequences 
for trapping ions in the dipole trap. We also describe 
the procedure for compensating stray electric fields, such 
that the resulting forces are smaller than the maximal 
forces in the dipole trap. Likely reasons for these stray 
electric fields are presented. Afterwards, the ramp-down 
procedure of the RF potential for the transfer of the ion 
into the dipole trap is described and discussed with the 
help of the stability diagram of a RF trap. Again, the 
influence of the static electric potential is highlighted. 
In the following, we describe the dipole trap setup. The 
method to adjust the polarization of the dipole trap beam 
to a purity of better than 1 : 1000 (intensity) is ex- 
plained, which is of importance to obtain the maximal 
potential depth of the dipole trap and minimal heating. 
Furthermore, a brief discussion of the heating effects in 
dipole traps is presented revealing that recoil heating is 
expected to be dominant in our dipole trap beam. 

In the following the experimental results are analysed. 
The analytic model used in Ref. 1 is derived. It assumes 
several simplifications, but is in agreement with the the- 
oretical expectation of recoil heating as the dominant 
heating effect. One of the simplifications is the cutoff 
of the trapping probability distribution in energy space, 
while in the experiment the cutoff happens in the posi- 
tion or potential energy space, respectively. An improved 
analytic model is introduced to overcome this simplifica- 
tion. Furthermore, the improved model allows for an 
estimate of the recoil heating rate as a function of tem- 
perature, which was previously assumed to be constant. 
Still, both analytic models simplify the dipole trap poten- 
tial as a spherically symmetric harmonic potential. We 
analyse the asymmetric shape of the potential consisting 
of dipole trap and static electric potential (including its 
defocussing effect) and quantify its depth as a function 
of the parameters of the static electric potential. Our 
setup offers only limited possibilities for measuring the 
parameters of the static electric potential, but according 
to our analysis these partly unknown parameters have a 
substantial effect on the overall trapping potential. In 
the best case scenario they effectively reduce the depth 
of the dipole trap by already 25 %. In the worst case 
scenario, the overall potential does not have a local min- 
imum, meaning that trapping becomes impossible. To 
refine the analysis, a Monte-Carlo simulation has been 
developed, which accounts for the asymmetry of the to- 
tal potential. Besides several other experimental details, 
the Monte-Carlo simulation also accounts for the Dopplcr 
cooling process in the RF trap and the transfer of the ion 
into the dipole trap. It is in good agreement with recoil 
heating as the dominant heating effect and supports the 
results of the analytic models for the experimental pa- 
rameters. However, it also reveals a limited validity of 
the analytic models for significantly different parameters. 
Finally, we conclude with a discussion of the experimen- 
tal limitations of our setup, propose improvements, and 
give an outlook of possible applications for optical trap- 



ping of ions. 

The article is organized as follows: In Sec. I we discuss 
the experimental procedures, their peculiarities, and dif- 
ficulties. In Sec. II, we provide a comparison of the dif- 
ferent models used to analyse the experimental results. 
Finally, in Sec. Ill, we present concluding remarks. 



I. EXPERIMENTAL PROCEDURE 

A single 24 Mg + ion is created by photo-ionization from 
a thermal atomic beam and Doppler cooled to a few mK 
(Doppler cooling limit: 1 mK) in the segmented linear 
RF trap [9]. We currently use the zone consisting of 
the three DC segments with the greater electrode-ion 
distance of Rq = 0.8 mm (Fig. 1). The radial frequencies 
in the RF trap result from the RF potential and are set to 
uj x .y « 2ir x 900 kHz corresponding to a potential depth 
°f U x ^ y R3 x 10 4 K, where denotes the Boltzmann 
constant. 



A. Axial Confinement by DC Electric Fields 

The axial frequency is tuned to uj z > 2ir x 100 kHz dur- 
ing the loading process of the RF trap by applying typi- 
cally positive DC voltages Vbc,out to the outer segments. 
The ion has a distance of approximately Zq = 1.5 mm 
to the edges of the outer DC segments. In principle, we 
would expect a focussing effect along the z-axis over the 
complete range of ±Zq due to the DC fields and a well 
depth on the order of U z ~tex 10 3 K. 

However, the axial potential of the RF trap shows 
strong anharmonicities, in particular, at a low axial fre- 
quency u) z R3 2it x 45 kHz, which is chosen during the 
experimental procedures (i.e., compensation of stray elec- 
tric fields and trapping in the dipole trap) described later. 
The images in Fig. 2 show ion chains at different axial 
frequencies. The symmetry for the high axial frequency 
indicates that a harmonic confinement is a good approx- 
imation for the respective parameters (Fig. 2a), while at 
the low axial frequency the shift of the center of mass 
and the different mutual distances between neighbouring 
ions visualize the axial anharmonicities (Fig. 2b). Due to 
the anharmonicities it is not possible to trap more than 
three ions at the low axial frequency chosen for Fig. 2b. 

The potential is characterized quantitatively in the fol- 
lowing: We extract the (absolute) positions of one, two, 
and three ions from Fig. 2b and measure the oscillation 
frequency of the single ion by resonantly exciting its ax- 
ial motion with a modulation of the voltage of a DC seg- 
ment. (We do not consider the oscillation frequencies of 
multiple ions, because the anharmonic potential inhibits 
sharp resonances then.) The six ion positions and the 
oscillation frequency are used to fit the parameters u; cu b. 
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Figure 2. Fluorescence images of one, two, and three ions 
in the RF trap. The radial frequencies are set to u) x , y « 
2ir x 870 kHz, while the axial trap frequency of a single ion 
amounts to (a) u z = 2n x (114.0 ± 0.5) kHz and (b) cu z = 
2tt x (40 ± 2) kHz, respectively. The dashed lines indicate the 
position of a single ion. The objective lens used to image the 
ions is intentionally defocussed to gain additional information 
about the ions' positions perpendicular to the imaging plane: 
The ion images will become larger (smaller), if the ion moves 
towards (away from) the CCD camera. Figure (b) reveals 
anharmonicities of the axial potential. 
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Here, L f DC,cub(^) denotes the cubic approximation of the 
the static electric potential, -Fboppicr the z-component of 
the (estimated) force corresponding to the light pressure 
of the Doppler cooling beams, and m the mass of the ion. 

The result of the best fit of Uuc,cub(z) is presented in 
Fig. 3. The best fit values for w cu b and z\ reproduce well 
the oscillation frequency of a single ion and its position, 
respectively. The best fit value z-z — — 135 um determines 
the distance of the local maximum to the local minimum 
and can be interpreted as a measure for the strength of 
the anharmonicities. The potential has a well depth of 
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Figure 3. Axial potential of the RF trap in cubic approxi- 
mation f/DC,cub and harmonic approximation [/DC,harm- The 
cubic approximation is derived from the positions of one, two, 
and three ions and the frequency of the axial oscillation of a 
single ion (Fig. 2b). The light pressure of the Doppler cool- 
ing beam is substracted. The axial potential has a maximum 
at a distance of 135 pm from the minimum and a depth of 
approximately fee x 0.5 K. 



The asymmetric potential is probably due to the load- 
ing process that coats the trap electrodes with Mg lead- 
ing to contact potentials (— 1.44 V for Mg on Au). This 
assumption is supported by the fact that a negative volt- 
age of —1.47 V has to be applied to one pair of outer DC 
segments for axial frequencies around u> z ~ 2ir x 45 kHz, 
while the voltage applied to the center segments is close 
to zero. This can only be explained by confining po- 
tentials on length scales smaller than Z$. (The ax- 
ial frequency due to the RF potential is estimated to 
~ 2tt x kHz and can be precluded as reason.) The appar- 
ent solution of shuttling the ions into a zone of the RF 
trap, where the electrodes are not coated by the atomic 
oven, is not an option in our setup, because the electrode- 
ion distance is smaller (Rq = 0.4 mm) and the dipole trap 
beam would produce significant stray light. 

The strong dependence of the axial potential depth on 
the axial frequency prevents a further reduction of the 
static electric potential below u z ~ 2ir x 40 kHz during 
the dipole trap experiments. In addition, the ions are 
prone to loss by collisions with components of the residual 
gas already at the axial well depth of k-Q x 0.5 K leading 
to a lifetime on the order of minutes. To increase their 
lifetime the axial frequency is increased to uj z > 2tt x 
100 kHz after the experimental procedures and during the 
loading process. 

Note that the axial potential can still be well approx- 
imated by a harmonic potential while trapping an ion 
in the dipole trap, because the displacement of the ion 
from the center of the RF trap remains small. In these 
experiments an anharmonicity of the total potential is 
due to the angle of 45° between the dipole trap beam 
and the z-principal axis of the static electric potential 
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(cmp. Sec. II B). 



B. Compensation of Stray Electric Fields 

In the presence of static electric fields the ion is dis- 
placed from the node of the RF potential. If the RF 
voltage and thus the stiffness of the RF potential is de- 
creased, the ion will be further displaced. Loading the ion 
into the dipole trap might fail, because the static electric 
fields can easily exceed the maximal forces of the dipole 
trap. 

Static electric fields have two major contributions: 
Firstly, the DC voltages applied to the DC segments lead 
to well controllable fields. Secondly, charges on dielectrics 
in the surrounding of the RF trap result in stray electric 
fields, which vary during the experiments. The optimiza- 
tion of the DC voltages and compensation of the stray 
electric fields is achieved with the ion as a sensor: The 
RF voltage is gradually ramped down and the ion's dis- 
placement is counteracted by appropriate DC voltages 
applied to the DC segments (see Fig. 1) and an addi- 
tional wire placed below the RF trap. For each itera- 
tion step of the compensation the RF voltage is ramped 
down only for a duration on the order of 0.5 s to main- 
tain the temperature of the RF circuit. Currently, we are 
able to perform the compensation without losses down 
to approximately 15 % of the initial amplitude of the 
RF voltage and thus to 15 % of the initial frequency 
(w Xiy w 0.15 x 2tt x 900 kHz = 2tt x 135 kHz). How- 
ever, trapping in the RF trap was not yet possible at 
< 10 % of the initial amplitude. The voltage applied to 
one of the inner DC segments is fine-tuned on the or- 
der of a few 100 uV. Due to the electrode-ion distance 
of Rq = 0.8 mm and geometric considerations this corre- 
sponds to a residual force on the order of 10 -20 N at the 
position of the ion. 

As stated above the RF potential of a finite linear RF 
trap has a small contribution to the axial confinement . A 
displacement of the ion during the ramp down of the RF 
voltage may also be observed in the axial direction. In 
this case, the axial position of the ion must be adjusted 
accordingly. However, this has to be done only once, 
because the node of the RF potential is fixed. 

We speculate that there are two main sources of stray 
electric fields in our setup, some of which require a com- 
pensation after each loading of an ion: Firstly, dielectrics 
in proximity of the ion are charged due to the photo- 
electric effect by stray laser light. The timescale for the 
discharge amounts to several minutes. Furthermore, as 
mentioned in Sec. I A, each loading process further coats 
the trap electrodes and thereby changes the contact po- 
tentials. The resulting changes of the axial potential 
(cmp. Fig. 3) are not relevant on timescales of weeks and 
the durations of the experiments presented in Sec. II, re- 
spectively. 

For quantum simulation experiments in this RF trap 
[10, 11] the fluorescence modulation technique from 



Ref. 12 has been used to optimize the DC fields and min- 
imize micromotion, respectively. However, after a micro- 
motion compensation with this technique a reduction of 
the RF voltage still displaces the ion. We speculate that 
thermal effects due to other operational parameters of 
the RF trap in the quantum simulation experiments lead 
to this discrepancy. That is why the fluorescence modu- 
lation technique might be adequate to minimize the mi- 
cromotion for stiff and constant confinements, but might 
be less suitable for a compensation of stray electric fields 
in the case of negligible RF confinements. 

C. Ramping Down the RF Voltage 

Besides of the compensation of stray electric fields at 
the position of the ion, the shape of the static electric 
potentials is important for a successful trapping in the 
dipole trap. The DC electric potential can be approxi- 
mated to quadratic order at the position r = (0, 0, 0) of 
the ion: 

U^-\mH,^+^ y ,y"+^). (2) 

It must satisfy Laplace's equation 

ZIUdcM = u%, + w 2 , + uj 2 z = 0. (3) 

Hence, for real us z > 0, at least one out of cu x / i y i must be 
imaginary resulting in a defocussing effect in the corre- 
sponding direction. Due to the restrictions in the reduc- 
tion of lu z discussed in Sec. I A, the defocussing is com- 
paratively strong in our current setup. The defocussing 
effect will not only overcome the weak confinement of 
the dipole trap, if the static electric potential is not ad- 
justed carefully. It also prevents trapping in the RF trap 
at small RF voltages, which are eventually applied while 
transferring the ion from the RF into the dipole trap. In 
the following we discuss this defocussing effect and the 
ramp-down procedure of the RF voltage with the help of 
the stability diagram of the RF trap. 

The radial motion of the ion in the RF trap (with- 
out considering an axial confinement at the moment) is 
described by Mathieu's equation [13, 14]. Whether trap- 
ping is possible, depends on the parameters q = jp^fp 

and a — ±a x t / y i — "T^^jp 2 "- Here, e denotes the ele- 
mentary charge, 4?x = 2V X 56 MHz the drive frequency, 
m s» 24 u the mass of the ion, Vrf the amplitude of the 
RF voltage, and Vbc,com a common DC voltage applied 
to all DC segments. The coordinates x' and y' corre- 
spond to a coordinate system rotated about the z-axis 
ifix> ly> = (fix T 6y)/v / 2 in this context; cmp. Fig. 1). The 
lowest order stability region for pairs (q, a) is shown in 
Fig. 4. Carelessly decreasing the RF voltage eventually 
leads to the loss of the ion before zero voltage is reached 
and is similar to the scan of a RF massfilter. 

However, it will still be impossible to approach arbi- 
trarily small RF voltages, while remaining in the sta- 
ble region, if the axial confinement and its defocussing 
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Figure 4. Lowest order stability region of Mathieu's equations 
of an infinitely long linear RF trap (without considering an 
axial confinement). The motion is stable for pairs (q, a) inside 
the gray shaded region. Outside this region the motion is 
only stable in x'-direction or j/'-direction or unstable in both 
radial directions. The arrow indicates a path in the stability 
diagram, when the RF voltage is decreased. It starts from 
a typical operational point of the RF trap. The motion of 
the ion will become unstable for any n ^ 0, before zero RF 
voltage is reached. 
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Figure 6. RF voltage of the RF trap during the experimental 
sequence. The amplitude is measured with a pick-up antenna 
inside the helical resonator, followed by filters to suppress 
extraneous RF signals from other sources in the laboratory. 
It is recorded by an oscilloscope, which is the reason for the 
noise. The final ramp-up of the RF voltage is slightly non- 
linear as the ramping speed does not meet the specifications 
of the RF synthesizer. The magnification (inset) shows the 
effect of the TTL-controlled RF-switch for times t < Topticai- 
The confinement due to the residual RF voltage during T op ti C ai 
can be neglected. 
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Figure 5. Example for lowest order stability region of Math- 
ieu's equations of a linear RF trap considering a specific de- 
focussing due to an axial confinement. We assumed some 
Voc.out 7^ and different geometry factors n x > 7^ K y i. It 
demonstrates that a decrease of the RF voltage will eventu- 
ally lead to an ion loss for any a-parameter (even a = 0) 
before zero RF voltage is reached. Additionally, it shows that 
Vrjc.out = need not be the optimal value to stay inside the 
stable region as long as possible. 



effect, respectively, are taken into account. It can be 
considered by introducing voltages Vbc,out 7^ Voccom 
at the outer DC segments. The a-parameters of Math- 
ieu's equation change according to a x i j y t — > a x i j y i — 

Kx'/y' 4e(VDC f j2 t jJf c,c ° m) , where K x >/ y > denote geometry 



factors [14]. Alternatively, the stable regions for the mo- 
tion in x'- and y'-direction, respectively, can be shifted 
for a specific Vbc,out, while the a-parameters remain un- 
changed (see Fig. 5). 

In general, the electrode geometry and the contact po- 
tentials due to the coating of the electrodes with Mg lead 
to a rotation of the x' and y' principal axes by an angle 
a referred to the coordinate axes and the geometry fac- 
tors k x i iyi are different. The values of lj x i and uj y i as 
well as the angle a are not experimentally accessible in 
the current setup and the contact potentials complicate 
predictions. A deeper analysis is presented in Sec. II C 
in the context of the Monte-Carlo simulation. We note 
that the discussion only considered an ideal quadrupole 
potential. Higher order contributions of the RF potential 
lead to instabilities within the theoretically stable regime 
[15], which can also result in the loss of the ion. 

In the dipole trap experiment the RF amplitude is 
ramped down linearly in time within 50 us for the trans- 
fer of the ion from the RF trap into the optical dipole 
trap and vice versa. In this way, the transfer is suffi- 
ciently slow with respect to the timescales set by any 
oscillation frequency. Still, instabilities are crossed com- 
paratively fast to mitigate their impact. At the same 
time the transfer is slow compared to the ring-down time 
of the helical resonator enhancing the RF voltage for the 
RF trap (1/e ring-down time: 0.5 us) and avoids over- 
shoots of the oscillating circuit. 

The ramp-down is achieved using the amplitude mod- 
ulation (AM) mode of the RF synthesizer driving the RF 
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Figure 7. Relevant part of the electronic energy level scheme 
of 24 Mg + . As 24 Mg + has no nuclear spin (7 = 0), there are 
only one ground level (S1/2) with angular momentum J = 1/2 
and two exited levels (P1/2 and P3/2) with angular momenta 
J = 1/2 and J = 3/2, respectively. The Mj sub- levels are 
shown. The dashed line indicates the cr + -polarized dipole trap 
beam with a detuning of A from the cycling transition. 



circuitry of the RF trap. The maximum attenuation by 
the AM amounts to approximately 99.75 %. An addi- 
tional TTL-controlled RF-switch allows to further atten- 
uate the RF voltage by a factor of 10 5 after ramping. 
Hence, the total attenuation of the RF voltage during 
the optical trap duration T opt ; ca i, i.e. the duration of 
maximal attenuation of the RF potential, is on the order 
of 10 8 . The confinement due to the residual RF poten- 
tial can therefore be neglected compared to the confine- 
ment of the dipole trap and the static electric fields. The 
ramping procedure is monitored on an oscilloscope using 
a pick-up antenna inside the helical resonator, followed 
by filters to suppress extraneous RF noise from other 
sources in the laboratory (see Fig. 6). 



D. Dipole Trap 

In the following we discuss the dipole trap into which 
the ion is transferred during the ramp-down of the RF 
voltage. The laser setup for the dipole trap beam con- 
sists of a 2 W fiber laser at 1118 nm with two consecutive 
second harmonic generation (SHG) stages [16]. The first 
SHG stage has an efficiency of around 70 % and the sec- 
ond SHG stage of around 40 %. This leads to a power ex- 
ceeding 500 mW at a wavelength of 280 nm with a highly 
astigmatic beam profile and a ratio of vertical to horizon- 
tal waist radii > 10 : 1. (At later times we even achieved 
a maximal power of around 600 mW.) The beam is 
shaped by a cylindrical and a spherical telescope. Further 
optical elements include an accousto-optical modulator 
to control its power, polarization optics, and a spheri- 
cal 1 : 10 telecope with pinhole to expand and clean the 



mode of the beam. The expanded beam is focussed by 
a spherical lens with a focal length of 200 mm onto the 
ion. Currently, we have up to P op ti C ai = 275 mW in the 
dipole trap beam. 

The laser can be detuned up to A = —2irx 300 GHz red 
of the resonance of the S1/2 <-> P3/2 transition of 24 Mg + 
with uj q = 2tt x 1 072.082 934 THz [17] and a natural line 
width of r = 2tt x 41.8 MHz [18] (see Fig. 7 for a level 
scheme). Thus, the maximal detuning corresponds to 
Z\w-7200r. 

During the dipole trap experiments a typical power of 
-Popticai = 190 mW is used in the dipole trap beam and 
a detuning of A ps —2tt x 275 GHz. The beam is fo- 
cussed on the ion and has a nearly Gaussian shape with 
a waist radius of u>o ss 7pm (1/e 2 radius of intensity). 
The polarization is tuned to a + (see Sec. IE). The first 
few off-resonantly scattered photons from the dipole trap 
beam pump the ion into the Sin(Mj = 1/2) level. As 
S 1/2 (M,7 = 1/2) «-> P 3/2 (M,7 = 3/2) is a cycling transi- 
tion, the ion can be treated as a two-level system. 

Considering solely the dipole trap potential we expect 
a depth of Uq = P °w^ ~ — X 38 mK and a max- 
imum force perpendicular to the beam of F ra d ~ 10 -19 N 
[19]. The corresponding trapping frequencies at the ap- 
proximately harmonic bottom of the dipole potential 
amount to w ra d = \J — AUq/ {mw^) s» 2tt x 165kHz in 
the radial directions. The confinement in the direction 
of beam propagation is Wf. ~ 2tt x 2 kHz and can be ne- 
glected. 

We note that, on the one hand, the beam profile of the 
dipole trap beam is not perfectly round due to experi- 
mental imperfections, but has waist radii of u>o,h « 5 pm 
in the horizontal direction (x-z-plane) and wo, v ~ 7 pm 
in the vertical y-direction. This is a relict of the beam 
shape from the second SHG stage and can hardly be com- 
pensated without sacrificing too much laser power. On 
the other hand, any deviation from a perfect Gaussian 
mode (due to, e.g. abberations, diffraction, . . . ) easily 
leads to a reduction of the peak intensity and thus a re- 
duction of the depth of the dipole trap potential. We 
think that the trap depth calculated for the round beam 
with wq = 7 pm with an measurement error (including 
both statistical and systematic errors) of 1 pm remains 
an appropriate estimate. 



E. Optimization of the Polarization 

Besides of the shape of the dipole trap beam, a 
pure polarization is important to allow for the maximal 
dipole trap depth and minimal heating (see the follow- 
ing Sec. IF). The polarization is adjusted to cr + referred 
to a quantization axis provided by a magnetic field of 
B = 5.5 G. The magnetic field is parallel to the propaga- 
tion direction of the dipole trap beam. The polarization 
is optimized using a 25 Mg + ion that has a nuclear spin of 
/ = 5/2. Hence its ground state is hyperfme-split into the 
S 1 / 2 (i 7 ' = 2) and Si/ 2 (F = 3) manifolds with an energy 
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Figure 8. Theoretically expected maximal total energy gain 
(d_E/dt)tot,max and maximal energy gain due to recoil heating 
(d_E/dt) s . max as a function of the detuning A for a Gaussian 
beam dipole trap (according to Ref. 20). The graphs are 
calculated for a waist radius wq = 7 pm, a potential depth 
Uo = —fee x 38 mK, an position in direction of the beam 
propagation zu = of the ion, and the radial position r with 
maximal heating rate. The gray shaded area is magnified 
in the inset, which allows to extract the energy gain for the 
experimental parameters. The two curves are almost identical 
such that for the given parameters recoil heating is expected 
to be dominant. 



difference of 2irH x 1.8 GHz. Neigbouring Mp sub- levels 
are Zeeman-shifted by ps 2tt x 3 MHz by the magnetic 
field. The dipole trap beam at a low intensity is tuned 
close to the S 1/2 (F = 3, Mp = 3) -H- P 3/2 (F = 4,M F = 
4) cycling transition. Polarization components other 
than tr + can lead the ion outside the cycling transition 
into the Si/ 2 (F — 2) ground state, which is off-resonant. 
Hence, the ion appears "dark" in state Si/ 2 (F = 2) and 
fluoresces "bright" in state Si/ 2 (F = 3). The optimiza- 
tion of the polarization is performed by analysing the 
fluorescence of the ion and adjusting polarization optics 
accordingly. Small coils producing magnetic fields per- 
pendicular to the main field are used align the quantiza- 
tion axis with the direction of beam propagation. This 
yields a a + purity on the order of typically 10 5 : 1 (inten- 
sity). In the experiments with 24 Mg + the polarization is 
not re-adjusted on short timescales, still allowing for a 
modest estimate for the purity of better than 1000 : 1. 

The Zeeman shift of the Mj-levels in 24 Mg + amounts 
to approximately ±2ir x 9 MHz, which is negligible com- 
pared to the AC Stark shift due to the dipole trap beam 
on the order of 2tt x 800 MHz. 



F. Heating Effects in Dipole Traps 

The analysis of the experimental data requires the 
knowledge about the relevant heating effects in dipole 
traps. Heating effects in dipole traps have been theoreti- 



Figure 9. Theoretically expected maximal total energy gain 
(d£'/dt)tot,max and maximal energy gain due to recoil heating 
(d_B/dt) s . max as a function of the detuning A for a standing- 
wave dipole trap created by retro-reflecting the Gaussian 
beam (according to Ref. 20). The graphs are calculated for a 
waist radius wo = 7 pm, a potential depth Uo = — x 38 mK, 
a radial position r = of the ion, and the position in direc- 
tion of the beam propagation Zk with maximal heating rate. 
The power P op ticai(^) of the dipole trap beam corresponds 
to a quarter of the power of Fig. 8 to account for the retro- 
reflection and interference in the standing wave. The gray 
shaded area is magnified in the inset. The total heating rate 
approximately doubles at the detunings used in the single- 
beam dipole trap experiments. 



cally examined in the 1980s in Refs. 20 and 21 for a two- 
level system with negligible velocity. The main heating 
effects were determined to be recoil heating and heating 
due to dipole force fluctuations. They are quantified in 
terms of a momentum diffusion constant 2D p , which re- 
sults in an average energy gain of dE/dt = D p /m for the 
trapped particle. 

The recoil heating is caused by a momentum diffu- 
sion due to absorption of laser photons and the subse- 
quent spontaneous emissions resulting in a random walk 
in momentum space. Each of the absorptions and spon- 
taneous emissions increases the expectation value of the 
variance of the momentum a 2 by (hk) 2 , where k — luq/c 
denotes the wavenumber of the laser photons. Thus, the 
energy gain can be written as (dE/dt) s — T s x 2E rec 
with the scattering rate r a and recoil energy E rcc — 
(hk) 2 /(2m) ps k-Q x 5uK. The recoil heating rate has 
its maximum F S)max = (Uo/h)(r/A) at the center of 
the dipole potential, where the intensity has its max- 
imum. For a potential depth of Uq — — k^ x 38 mK, 
a detuning of A — ~2tt x 275 GHz, and a waist ra- 
dius of wo = 7um we obtain -T Sjmax ps 750 ms -1 and 
{dE/dt) a>max w fc B x 7.7mK/ms. 

Dipole force fluctuations are due to the opposite signs 
of the dipole forces for an ion in the states S1/2 (Mj — 
1/2) and P3/2 (Mj = 3/2). In the dressed-state approach 
of Ref. 21 spontaneous emissions between levels subject 
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to different dipole forces lead to a momentum diffusion 
and energy gain, respectively. The diffusion constant 
is proportional to the square of the intensity gradient. 
Therefore, the energy gain due to dipole force fluctua- 
tions is highest close to the gradients of the dipole po- 
tential. 

The full mathematical treatment for a two-level atom 
at rest can be found in Ref. 20. In Fig. 8 the total en- 
ergy gain (di£/di)tot as a function of the detuning A is 
compared to the energy gain solely due to recoil heating 
(dE/dt) s . The potential depth Uq is kept constant by 
adapting the power P op ticai(^) as a function of the de- 
tuning. The difference (dE/dt) to t — (dE/dt) s amounts 
to only k B x 0.02 mK/ms at A = -2?r x 275 GHz. Thus, 
at the typical experimental parameters recoil heating is 
expected to be dominant. 

However, for a standing wave dipole trap created by 
retro-reflecting the Gaussian beam, the heating due to 
dipole force fluctuations has to be considered (see Fig. 9) . 
This is due to the comparatively large intensity gra- 
dients in direction of the standing wave (the potential 
now has a modulation on length scales A/2 -C wq). 
The power P op ticai(^) was reduced by a factor of 4 
compared to Fig. 8 to maintain the potential depth of 
Uq = —fee x 38 mK. The maximal total energy gain in- 
creases by roughly a factor of two at our experimental 
parameters compared to the trap consisting of a single 
Gaussian beam: (d£ , /dt) to t.max ~ &b x 15 mK/ms. 

The dependence of the heating rate due to dipole force 
fluctuation on the Rabi frequency /2r can be estimated 
from Eq. (4.6) of Ref. 2f . ft yields D p <x f2^ for detunings 
A ^> As Uq oc J?j|, we expect an increase of the heat- 
ing rate by a factor of f 6 for a doubling of the potential 
depth at A = —2tt x 275 GHz. The exact value according 
to Ref. 20 is (d-E/di)tot,max ~ &b x 220 mK/ms, which 
emphasizes the strong dependence on the intensity. This 
adds difficulty to dipole trap experiments with a standing 
wave. 

Further heating in the dipole trap can be caused by 
an impure polarization. The coupling of the dipole trap 
beam with a + polarization to the Si/2(Mj = —1/2) -H> 
P3/2 (Mj = 1/2) transition is a factor of l/\/3 weaker 
than to S1/2 (Mj = 1/2) o P 3 / 2 (-Mj = 3/2). Polariza- 
tion components other than a + can scatter the ion into 
the Si/2(M,/ = —1/2) level. A fluctuation between the 
S 1/2 (Mj = 1/2) and S 1/2 (Mj = -1/2) levels leads to 
force fluctuations and thus a momentum diffusion (see, 
e.g., [22]), which is similar to the diffusion caused by 
dipole force fluctuations above. Additionally, the poten- 
tial depth in the Si/ 2 (Mj — —1/2) level is reduced to 
Uq /3, which can lead to an earlier loss of the ion. How- 
ever, we expect this effect to be negligible in our experi- 
ment (see Sec. IE). Heating due to intensity noise in the 
dipole trap beam [23] is also estimated to be orders of 
magnitude smaller than the recoil heating rate. 



II. EXPERIMENTAL RESULTS AND ANALYSIS 

The aim of the measurements is to obtain an estimate 
for the trapping probability P of an ion in the dipole 
trap for different parameters T op ti C ai and -Popticai- P is 
determined as the mean number of successful trapping 
attempts with approximately 40 ions for each set of pa- 
rameters. For example, this corresponds to more than 
200 trapping attempts according to the above scheme for 
^optical = 1ms and P op ticai = 190 mW (cmp. Fig. 10). 

Results of the dependence of the trapping probability 
on Topical and Popticai have been published in Ref. 1. 
To summarize, the half-life in the dipole trap can be es- 
timated to ?op^? a ] = 2.4 ms and to be mainly limited 
by the intrinsic recoil heating by the dipole trap beam. 
Furthermore, a measurement with P op ticai = mW re- 
sulting in zero trapping probability confirms that there 
is no other relevant three-dimensional trapping potential 
during T optical . 

A. Analytic Models 

The experimental data presented in Ref. 1 have been 
analysed using a simple analytic model. In this section, 
the analytic model (cmp. Ref. 1) is derived, and the im- 
pact of the simpilifications is discussed. Based on the 
conclusions an improved analytic model is developed. 

For the first analytic model the total potential con- 
sisting of dipole trap and static electric potential is ap- 
proximated as a spherically symmetric harmonic poten- 
tial with frequency Wharm^ 

U h3 , Tm (f) =U + ^mw 1 2 larm r 2 . (4) 

(Note that Uq < by definition.) In the model the parti- 
cle energies E are subject to the Boltzmann distribution 
of this potential: 

f(E) = ^(E-U^e-^-^ with = -L (5) 
2 fc B J 

where T denotes the temperature. 

A particle will be considered as trapped, if it has a 
total energy E < after the dipole trap duration T opt i ca i. 
That is why the model is referred to as the energy- cutoff 
model in the following. The trapping probability for a 
particle is obtained as [24] 

Pcn(0 = £ dE f(E) = 1 - 0* + £ + e"€ (6) 

with £ := P\U \. 

To model the dipole trap experiments, the temperature 
is assumed to increase linearly in time t with heating 
rate R := dT/dt from an initial temperature To: T(t) = 
To + Rt. Both the potential depth \Uo\ ='■ \uq | Popticai 
and the heating rate R =: i? Popticai (see Sec. IF) are 
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Figure 10. Comparison of the energy-cutoff and radial-cutoff 
model: The optical trapping probabilities P en according to 
the energy-cutoff model and P ra d according to the radial- 
cutoff model are shown as functions of the dipole trap du- 
ration Topticai- The points represent the measured optical 
trapping probabilities and the error bars the standard devi- 
ation of the mean [1]. The parameters for which the curves 
of the analytic models are shown originate from a joint fit 
to the experimental data of this figure and Fig. 11. While 
the curves are almost identical, the best fit values for the ini- 
tial temperature and the heating rate are approximately 20 % 
greater for the radial-cutoff model. Experimental parameters: 
wo = (7± l)um, A = -2tt x 275GHz, P optica i = 190mW, 
and Lu z = 2n x 47 kHz. 



assumed to be proportional to the power of the dipole 
trap beam. Hence, we obtain 



I WO I ^optical 



k B (T + R Popticai t) 



(7) 



In the experiments the heating in the direction of beam 
propagation is larger, because the absorption of photons 
always occurs in this direction as opposed to the emission 
of photons. However, the total potential in the dipole 
trap experiments is anharmonic (see also Sees. II B and 
II C), such that we expect cross-dimensional mixing. This 
justifies the assumption of an isotropic heating in the 
model. 

Next, we substitute £ in Eq. (6) with Eq. (7) and use 
the predicted value \uo\ = &b x 38mK/190mW. A joint 
fit of P en for the data of Figs. 10 and 11 yields parame- 
ters T = 4.6 mK and R = 3.9 mK/ms x P op ticai/190 mW. 
Here, the uncertainty of the waist of ±lum generates an 
error of |ito| and, hence, of the fitted parameters of 30 %. 
The time-derivative of the relation E = 3/cbT + C/o which 
holds for a three-dimensional harmonic potential yields 
dE/dt = 3/sbP- Hence, the value of the heating rate from 
the fit yields an energy gain dE/dt = fee x 11.7mK/ms, 
which is a factor of 1.5 higher than the estimated re- 
coil heating rate at the center of the dipole trap (cmp. 
Sec. IF) [25]. The following discussion will show that 
the experimental results are still in agreement with re- 
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Figure 11. Comparison of the energy-cutoff and radial-cutoff 
model (analogously to Fig. 10): The optical trapping prob- 
abilities P en and P ra d according to the models are shown as 
functions of the power Popticai of the dipole trap beam. Exper- 
imental parameters: wo = (6. 5± 1.0) urn, A = — 2-7rx300GHz, 
Topticai = 0.5 ms, and lo z = 2-zr x 41 kHz. 



coil heating as the relevant heating effect predicted in 
Sec. IF. 

Now we further analyse the energy-cutoff model. A 
major simplification is the cutoff of the probability dis- 
tribution in energy space. However, in the experiment 
the cutoff happens in position space or equivalently in 
potential energy space. Additional energy of up to \Uq\ 
can be stored in the angular momentum of a trapped 
particle. The energy-cutoff model is extended to a radial- 
cutoff model taking this into account by following a cal- 
culation in Ref. 26. The probability for a particle to be 
trapped at time t is defined as the probability to have 
an orbit with potential energy below zero. This is still 
an approximation, because a small fraction of the parti- 
cles is on orbits with potential energy larger than zero at 
intermediate times, but evolves back to lower orbits at 
a later time. These particles are not treated as lost in 
the radial-cutoff model. The possibility to evolve back to 
lower orbits arises from the diffusive nature of the heating 
process. 

In spherical coordinates the energy of the ion is given 

by 



E 



1 

2m 



I? 



+ U hablm (f), 



(8) 



where p r — mdr/dt denotes the momentum in radial 
direction and L the angular momentum. The maximal 
radius of an orbit is obtained for p r = 0: 



(E - U ) + y/(E-U ) 2 -(oJ haxm Ly 



1/2 



(9) 
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We define a characteristic function 



0{E,L) 



1 if £4armO m ax) < 

else 



(10) 



which is 1 only if the potential energy of an orbit does not 
exceed zero. For a harmonic potential the characteristic 
function can be rewritten as 

1 if U < E < 



G(E, L) = { 1 if < E < -U a and L > ^ J °f 2 ■ 
else 

(11) 

Hence, if a particle has an energy Uo < E < 0, it will 
always be trapped in the potential. If it has an energy 
< E < — Ua, it may (but need not) be trapped, because 
in addition to the energy of up to \Uq\ that can be stored 
as potential/radial kinetic energy, there is also kinetic 
energy of up to |Z7o| related to the angular momentum. 

With some effort, the trapping probability can be cal- 
culated analytically as a phase space integral: 

iW0 = / dnW(E)0(E,L) = l-e- 2 S-2Ze-S (12) 

with the infinitesimal phase space volume dJ? = d rd p 
and the density of states per phase space volume 



W(E) 



( 



2tt 



-P(E-U ) 



(13) 



The comparison of the trapping probabilities of the two 
analytical models shows P r ad(£) ~ Pen(l-2£). 

Again we substitute £ in Eq. (12) with Eq. (7) and 
fit -Prad to the measured data analogously to the fit of 
the energy-cutoff model. The best fit results for both the 
heating rate as well as the initial temperature are approx- 
imately 20 % higher than for the energy-cutoff model, 
which could be expected due to P ra d(£) ~ Pen(l-2£)- 
This shows that energy stored as angular momentum is 
of limited relevance. The fitted curves corresponding to 
both models are compared in Figs. 10 and 11 and are 
almost identical. 

To understand the deviations between experimental re- 
sults and analytic models, we have to further investigate 
the simplifications in the latters: Both models simplify 
the total potential of the trap as an isotropic harmonic 
potential. In contrast, the total potential in the dipole 
trap experiments is asymmetric and its well depth is de- 
creased in some directions due to the defocussing effect 
of the static electric potential. However, both analytic 
models still agree with the experimental results within a 
factor of two and are useful for rough estimates of ex- 
pected trapping probabilities. The potential shape and 
well depth will be further discussed in Sec. II B and the 
Monte-Carlo simulation presented in Sec. II C will finally 
overcome most simplifications of the analytic models. 

Before we conclude this section, we extended the calcu- 
lations of the radial-cuttoff model to analyse the scatter- 
ing rate. In the dipole trap potential it is proportional 



to the intensity and thus the potential itself: P s (r) oc 
— U(f). For the harmonic potential in Eq. (4), the corre- 
sponding scattering rate becomes 



Ps,harm if) 



.max arm 

(r)/U for £/harm(r) < 







else 



(14) 

The average scattering rate can be calculated analogously 
to P rad (0 : 



-^s, harm (0 



l 



Prad(e) 



df2W(E)0(E,L)r s 

.harm (r) 



(15) 



-r (^ 3-(2g + 3)c- 2 €-2g(g + 2)e-€ \ 

s ' max ^ 2£(1 - e- 2 ? - 2£e-* ) )' 1 ' 

For zero temperature this yields P s ,harm = Ps,max- I n 
the limit k^T oo the avarage scattering rate decreases 

to p. 



s.harm 



-> P s , max /2. At k B T = \U \/3 (correspond- 



ing to E — 0), it reaches P Si ha 



OMR 



This 



result shows that the avarage scattering rate is always 
greater than P s , max/2 (as opposed to the rough estimate 
in Ref. 1). A simplification in this calculation still is the 
isotropic decrease of the scattering rate in all three di- 
mensions, while in the experiment it remains constant 
along the propagation direction of the dipole trap beam 
to a good approximation. 



B. Analysis of Potential Shape and Well Depth 

The potential of the dipole trap close to the waist of 
the laser beam (distance to waist small compared to the 
Rayleigh length zq) can be approximated by 



^optical (r) ~ C^oexp 



(x — zy 



y 



(17) 



where the dipole trap beam propagates in direction 
k = (1,0, l)/v2- The trapping frequency at the ap- 
proximately harmonic bottom of this potential is CLi ra( j = 
^J—AUq/ (ttt-Wq). To simplify matters, we assume that 
the principal axes of the static electric potential £/dc co- 
incide with the coordinate axes (i.e., the angle between 
principal axes and coordinate axes is a = ° ) : 

Uuc(f f ) = ^m(u 2 x x 2 +cjy +u 2 z z 2 ). (18) 

The total potential is the superposition of these two po- 
tentials: U(r) = [/optical (r) + Ur>c(r)- 

The well depth of U{r) differs for different directions, 
because Uy>c{t) has a defocussing effect perpendicular to 
the z-axis. At some points the potential U(r) turns from 
focussing to defocussing for values U(r) < 0, i.e., the well 
depth is reduced compared to the depth \Uq\ of the dipole 
trap potential only. We are interested in the minimal 
particle energy U mpe which is neccessary to escape from 
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Figure 12. Equipotential surfaces U(r) — ?7 mpeima x of the 
total potential. The parameters of the graphs correspond 
to the experimental parameters Uo = — &b X 38 mK, ui z = 
2ir x 47 kHz, and wo = 7um. The frequency in {/-direction 
is (a) u? y = c*4 imax , (b) u* = (1.2cj y , max ) 2 , and (c) <j? y = 
(0.95 Lu y ,max) 2 , respectively (see text). The points defining 
(7 ra pe,max in (a) appear as "cut" trough the egg-like surface 
enclosing the trap centre (indicated by the arrow). Note that 
Wj mal is negative, i.e., an increase of the absolute value of uiy 
leads to a stronger defocussing effect in {/-direction (b), while 
a decrease leads to a stronger defocussing in a direction in the 
z-z-plane (c). 



the trap. (The minimal well depth can then be written as 
\Uq — t/ mpo |.) Any particle with an energy E < U mpe will 
be trapped independently of its trajectory. A particle 
with larger energy may be trapped for a certain duration 
depending on its trajectory. 

To calculate U mpc we derive the critical points r cr it 
of the potential, which fulfill VJ7(r)|f crit = 0. By sym- 
metry there are at least two such points ±r mpc with 
C/(±r mpe ) = U mpc in addition to the local minimum at 
r*o = (0,0,0) with U(fo) = Uq (as long as the poten- 
tial provides a three-dimensional confinement). We give 
solutions for some special cases in the following. 

Except for u> x and uj y all parameters of the poten- 
tial are sufficiently well known in the experiments (cmp. 
Sec. I). Therefore, we analyse U mpe as a function of u) v in 
the following, while w x = [— (u)y + W 2 )] 1 / 2 is varied along 
with ujy to meet Laplace's equation. An interesting case 
is the maximum of U mpe , because it sets the deepest trap 
we can expect. It is obtained for w 2 max = (1 — v^u; 2 
and amounts to J7 mpe ,max = bUo, with 

b=p [l-ln(p)l and p= (19) 

W rad 

By inserting oj z = 2ir x 47 kHz and w ra d — 2n x 165 kHz 
for the experimental parameters (cmp. Fig. 10) we obtain 
a defocussing frequency in y-direction of w 2 max = — (2n x 
40kHz) 2 , b « 0.23, and [/ mpc , max = -k B x 8.6mK. This 
means that particles see a minimal well depth amounting 
to only \Uo — U mpe \ = k-g, x 29.4mK in some directions 
even for an "optimal" choice of u) y . The well depth is 
reduced by already a quarter compared to \Uq\ = fee x 
38mK. 

The potential shape for uj y = Wy jlnax is special insofar 
as there is an infinite number of critical points r mpe lying 
on a "ring". This is graphically illustrated in Fig. 12a, 
where the equipotential surface for U(r) = U mpe max 
shows these points as "cut" (indicated by an arrow) 
through the egg-like surface containing the trap centre. A 
particle with energy E — C/ mpe , m ax can only escape from 
the trap through this "cut". For a larger defocussing in 
y-direction we obtain U mpc < J7 mpe , m ax- Two "holes" in 
the equipotential surface U(r) = U mpe>max form on the 
?/-axis (see Fig. 12b). Analogously, a larger defocussing 
in x-direction leads to two "holes" in the x-z-plane (see 
Fig. 12c). In the cases of Figs. 12b and 12c, a particle 
with energy E < U mpe)max can escape from the trap in 
direction of the "holes". 

Another interesting case is the value of oj y , where the 
total potential ceases to provide a three-dimensional con- 
finement. Two cases have to be distinguished: Firstly, 
the defocussing in y-direction trivially becomes larger 
than the focussing effect of the dipole trap beam for 
ujy = — w 2 ad R3 — (2tt x 165 kHz) 2 (again for the exper- 
imental parameters; cmp. Fig. 10). Secondly, for the x- 
direction the total potential becomes repelling for a (non- 
trivial) value of w 2 w — (27rx 43.5 kHz) 2 corresponding to 
u}y sa — (2n x 17.8 kHz) 2 . The higher sensitivity to a defo- 
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Figure 13. Probability distributions f(E) of the particle 
energies E during Doppler cooling in the RF trap. The 
curves show the results from the Monte-Carlo simulation 
and related fits of the Boltzmann distribution of the three- 
dimensional harmonic potential. The frequency splitting is 



= 2-7T x 2.5 kHz and the tilting angle to 1 — 1 



set to oj x i — 

a = 35 ° . The simulation is performed for a cooling duration 
of 2 s with time steps of At = 0.1/-T. "(i)" denotes the param- 
eters used during Doppler cooling corresponding to Fig. 10: 
u z = 2tt x 47 kHz, u x > = 2ir x 950 kHz, P DL = 8.6 uW, and 
ftiLdct = 270 uW. The axial displacement of the ion due to 
the light pressure of the Doppler cooling beams amounts to 
Az — 8.8 um. The radial displacement is of the order of few 
nanometres and can be neglected. The energies E are referred 
to the harmonic oscillator shifted by Az. The curves "(ii)" cor- 
respond to Fig. 11: luz = 2ir x 42 kHz, oj x i = 2n x 880 kHz, 
Pol = 2.1 uW, and PoLdet = 138 uW. These parameters 
yield a displacement due to light pressure of Az = 2.8 pm. 
The curves of the corresponding Boltzmann distributions fit 
well to the simulated curves. The temperatures according to 
the fits of the Boltzmann distribution are T (i) = 3.1 mK for 
the parameters of (i) and T' a ' = 2.7 mK for the parameters 
of (ii). 
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Figure 14. Comparison of the probability distributions f(E) 
of the particle energies E during Doppler cooling in the RF 
trap for (a) different splittings uj x i — u y i and (b) angles a. 
The curves show the results from the Monte-Carlo simulation 
with the same parameters as used in Fig. 13 for parameters 
(ii). For small frequency splitting u) x i — ui y i or small angles a 
Doppler cooling becomes inefficient. 



cussing in the x-direction can be explained by the angle 
of 45 ° between the dipole trap beam and the £-axis. 

The discussion reveals that the absolute value of the 
potential depth assumed in the analytic models (see 
Sec. II A) is overestimated, while the scattering rate is 
still determined by the laser intensity in the dipole trap 
beam. However, a certain C/ mpc does not imply that par- 
ticles with larger energy are eventually lost from the trap 
within T op ticai- Therefore, the analysis of the impact of 
the anisotropy of the potential requires a more realistic 
model, which is given in the following Sec. II C. 



C. Monte-Carlo Simulation 

In order to overcome the limitations of the models from 
Sec. II A, a rate equation Monte-Carlo simulation has 
been developed. Therein, the ion has a classical position 
f and velocity v and a discrete state i out of the level 



scheme depicted in Fig. 7. The simulation consists of 
two parts: Firstly, the Doppler cooling process in the RF 
trap is simulated. This simulation does not only allow 
to analyse the impact of the static electric potential on 
the efficiency of the Doppler cooling process in the RF 
trap. It also provides the initial values (position, veloc- 
ity, internal state) of the ions for the simulation of the 
dipole trap experiments and thus no free parameter for 
the initial temperature is required. Secondly, the dipole 
trap experiments are simulated. Therein, the ions are 
transferred from the effective potential of the RF trap 
into the dipole trap, in which they are exposed to re- 
coil heating. Algorithms from Ref. 27 were used to in- 
tegrate the equation of motion (Runge-Kutta-Fehlberg 
(4, 5) method) and to generate pseudo-random numbers 
from different distributions. 

For the Doppler cooling simulation the potential con- 
sists of the superposition of the time-averaged RF poten- 
tial and the static electric potential. For both individual 
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potentials one principal axis coincides with the z-axis, 
while the other principal axes e x > and e y > can be rotated 
by an angle a referred to the x- and j/-axis, respectively. 
(In general, the rotation angle of the time-averaged RF 
potential will be different from that of the static electric 
potential. The rotation angle of the total potential will by 
somewhere inbetween. This is neglected in the following.) 
The relation ui^,, +wj)< +^ = must hold for the frequen- 
cies of the static electric potential. For the time-averaged 
RF potential the relation w^ RF ' > — w^ff^ — = holds, 

where w^f" 1 ^ denotes the largest of the three frequencies 
urf, and uA R fL the two smaller ones. A small axial 

x/y/z <,l/2 

confinement CJ*^ ~ 2tt x kHz is due to the geometry of 
the electrodes. The axial frequency of the total potential 
is still well described by uj z (the correction due to Wz 
can be neglected) and the radial frequencies are split by 
few 2tt x kHz. Unfortunately, both the exact splitting 
of the radial frequencies as well as the angle a are ex- 
perimentally difficult to access in the current setup. The 
reliability of numerical calculations of the potentials of 
the RF trap is restricted due to the difficulties discussed 
in Sec. I A. 

In the Doppler cooling simulation both superimposed 
Dopplcr cooling lasers DL (with detuning of Z\dl ~ 
~r/2 from the S1/2 O P3/2 transition to cool to a 
low temperature) and DLdet (ZiDLdot ~ — 11-T to re- 
cool occasionally heated ions) are included. They are 
focussed onto the ion under an angle of 22.5° to the z- 
axis (fc DL = (0,sin(22.5°),cos(22.5°)), see Fig. 1) and 
have waist radii around 90 um. The polarization of the 
laser beams is the same and we assume a 7r-polarization 
in the simulation, as it is the main component. 

Each time step of the simulation consists of 

1. a propagation of the ion for a duration At in the 
potential, 

2. an absorption or induced emission of a photon 
from/into each Doppler cooling beam with a cer- 
tain probability, and 

3. a spontaneous decay with isotropic emission of a 
photon (for excited states i) with a certain prob- 
ability. 

The results of different Doppler cooling simulations are 
shown in Fig. 13. The powers of the Doppler cooling 
lasers -PoL/DLdct corresponding to Fig. 10 are larger (on- 
resonance saturation parameter for the laser DL so ~ 0.3) 
than the powers corresponding to Fig. 11 and Doppler 
cooling becomes slightly less efficient. The light pressure 
excerts a mean force of ~ 3 x 10~ 20 N onto the ion and 
leads to a significant displacement Az = 8.8 um. The re- 
sulting probability distributions deviate only lightly from 
the Boltzmann distributions of the three-dimensional 
harmonic oscillator, which are shown for comparison. 

In Fig. 14 the results of Doppler cooling simulations 
are compared for different splittings uv — uj y i and angles 
a, respectively. Doppler cooling becomes more efficient 
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Figure 15. Monte-Carlo simulations of the dipole trap ex- 
periments (cmp. Figs. 10 and 11) for an angle a = 35 and 
different frequencies uiyi. The initial values for the position, 
velocity and internal state are obtained from the respective 
results of the Doppler cooling simulations depicted in Fig. 13. 
The points represent the measured optical trapping probabil- 
ities and the error bars the standard deviation of the mean. 
The Monte-Carlo simulations agrees with the measured trap- 
ping probabilities for the depicted parameters. 



for large angles a and large splittings of the radial fre- 
quencies. While the impact of the angle is large, different 
splittings lead to only a small change as long as they are 
larger than ~ 27rx 1 kHz. Hence, the result of the Doppler 
cooling is comparatively insensitive to any (plausible) es- 
timates of the splitting of the radial frequencies. 

These dependencies of the Doppler cooling results are 
consistent with theory: Reducing the power of a Doppler 
cooling beam, reduces the mean energy of the particle, 
as long as external heating effects remain negligible. As 
the cooling beams are sent in along one direction only, 
they must not be (close to) perpendicular to any of the 
principal axes of the trap potential for Doppler cooling 
to be effective in all dimensions. This also implies that 
the frequencies of the total potential must be pairwise 
different. 

Next, the dipole trap experiments are simulated. The 
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initial conditions (position, velocity, and internal state) 
are picked from results of the Doppler cooling simula- 
tion. Initially, the potential consists of the one used dur- 
ing Doppler cooling superimposed with the dipole trap 
potential L/ O pticai(^0 &s in Sec. II B, but considering the 
slight asymmetry in the waist radii in horizontal and ver- 
tical direction (wq h = 0.7 Wo v = 0.7 Wq). In the exper- 
iment the dipole trap beam is focussed on the ion by 
observing a decrease of the Doppler cooling fluorescence 
due to the AC Stark shift caused by the dipole trap beam. 
The power of the dipole trap beam and thus the decrease 
of the fluorescence is kept as small as possible, such that 
the dipole trap beam is effectively focussed onto the ion 
displaced by the light pressure caused by the Doppler 
cooling beams. That is why the dipole trap potential in 
the simulation is also displaced by Az from the center of 
the potentials of the RF trap. The frequencies <^^fyi/ z of 
the time-averaged RF potential are linearly ramped down 
to zero within the first 50 us of the simulation. Both w y < 
(or u x i) and a are free parameters as in the Doppler 
cooling simulations. 

The trapped ion randomly absorbs photons of the 
dipole trap laser with the scattering rate r s (r). As 
r s (r) <C r, every absorption is followed by an imme- 
diate, isotropic spontaneous decay. The ion will be con- 
sidered as lost from the trap, if its distance to the origin 
exceeds several waist radii wq. Hence, we can be sure 
that it left the trapping region of the dipole trap and 
within roughly 10 us <C Topical it will also be accelerated 
out of the region, where the RF trap could re-capture 
it. The trapping probability P for different parameters 
^optical and Popticai is determined as the fraction of ions 
that were not lost. 

In Fig. 15, the results are illustrated for an angle 
a = 35 ° and different frequencies uj y i . The best agree- 
ment with the measured optical trapping probabilities is 
obtained for w 2 , k, ~(2tt x 50kHz) 2 (and similarly for 
cjy, R3 — (2ir x 35 kHz) 2 , which is not shown). For both 
simulations the maximal trapping probability is reached 
for ujy, ~ — (2n x 40kHz) 2 . This is consistent with the 
expectations from Sec. II B. (Note, however, that the po- 
tential has been simplified in these analytic calculations 
and they strictly hold for a — 0° only.) The simula- 
tion of the lifetime measurement (Fig. 15a) shows that 
the comparatively large displacement Az > w of the 
ion does not prevent it from being trapped successfully. 
The comparatively slow ramp-down of the RF potential 
allows the ion to be displaced mainly along the propa- 
gation direction of the dipole trap beam and loose most 
of the potential energy related to the displacement Az 
from the center of the static electric potential. The sim- 
ulation of the optical trapping probability P versus power 
of the dipole trap beam P op ti C ai (Fig. 15b) shows a com- 
paratively sharp increase for small powers and P almost 
reaches 1. The deviation from the measured probabili- 
ties can most likely be explained by beam pointing in- 
stabilities in the experiment, that smoothen the sharp 
increase and reduce the maximal trapping probability. 
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Figure 16. Monte-Carlo simulations of one dipole trap exper- 
iment (cmp. Fig. 10) for a frequency uiyi = —(2iv x 50kHz) 2 
and different angles a between principal axes of the static 
electric potential and the axes of the coordinate system. The 
initial values for the position, velocity and internal state are 
obtained from the same Doppler cooling simulation of Fig. 13 
for all three curves. Note that this allows for a better com- 
parison of the influence of a on the dipole trap experiment, 
but assumes a better Doppler cooling result for small angles 
than predicted by the simulations in Fig. 14. The trapping 
probability shows a strong dependence on the angle a. 



These beam pointing instabilies are mainly caused by 
the convection near the warm magnetic field coils of the 
quantization field. 

Within a certain range of parameters both the Doppler 
cooling and the dipole trap experiments are quite insen- 
sitive to the actual values of the free parameters. The 
lifetime experiment, for example, is consistent within 
the error bounds of the experimental parameters with 
iu 2 y , = -(2tt x [30 kHz... 60 kHz]) 2 at an angle a = 35°. 
However, for frequencies w 2 , outside this range the trap- 
ping probability is significantly reduced or even zero. The 
optical trapping probability P drops for some specific an- 
gles a towards zero (see Fig. 16 for a = 45°) or is en- 
hanced (see Fig. 16 for a — 20°), as long as the same 
Doppler cooling efhcieny is assumed for all simulations. 

The Monte-Carlo simulation still does not account for 
some effects that most probably also cause a decrease of 
the optical trapping probability, e.g. micromotion during 
the Doppler cooling and ramp-down of the RF potential, 
pointing instabilities of the dipole trap beam including an 
overall drift, residual stray electric fields, an arbitrarily 
rotated static electric potential (including the z-principal 
axis), or a reduction of the peak intensity of the dipole 
trap beam deviating from a perfect mode due to diffrac- 
tion, abberations, etc. Still, it overcomes the most se- 
vere limitations of the analytic models. Its results are in 
agreement with the measured trapping probabilities and 
support recoil heating as the relevant heating source in 
the dipole trap experiments. Additionally, the Monte- 
Carlo simulation allows to reveal the large impact of the 
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radial frequencies (w^, , uiy, ) of the static electric poten- 
tial and of the angle a of the principal axes. The inves- 
tigations also show that it is unlikely that the trapping 
probability can be significantly improved in the current 
setup, at least, without cooling of the ion in the dipole 
trap. The Monte-Carlo simulation is expected to also 
be a useful tool for future experiments to allow for both 
quantitative predictions and analyses of measurements. 

III. OUTLOOK 

Our experiment leaves room for various technical im- 
provements. First of all, the recoil heating rate must 
be reduced to allow for longer trapping durations. This 
can be accomplished by a dipole trap beam with larger 
detuning A and higher laser power. Laser systems as 
described in Ref. 28, a laser at 1120 nm [1, 29], or sys- 
tems at other prominent visible or infrared wavelengths 
(e.g. 532 nm or 1064 nm) are established candidates for 
this purpose and available as turn- key systems. These 
wavelengths provide further advantages compared to the 
ultraviolet wavelength used in our experiment: Optical 
fibres are available that allow for short beam paths in 
air and hence stable beam pointing and beam shapes. In 
contrast, the length of the beam path in the presented ex- 
periments amounts to several metres and the beam shape 
of the dipole trap suffers from the strong astigmatism of 
the second SHG. Stray electric fields drifting on short 
timescales due to charges produced by stray laser light 
via the photoelectric effect should be mitigated for long 
wavelengths. Obviously, a reduction of the distance of 
the windows of the vacuum chamber to the trap would 
ease to achieve small and robust waist radii. 

The RF trap has proven to be a useful tool to allow 
for the sensing and compensation of stray electric fields 
and loading an ion into the dipole trap, although in prin- 
ciple an ion could also be created by ionizing an already 
trapped neutral atom. The RF trap used in the presented 
experiments was originally optimized for quantum sim- 
ulation experiments [10, 11] requiring high secular fre- 
quencies. This results in small distances of the ion to the 
electrodes. As opposed to that, the dipole trap experi- 
ments would benefit from a trap optimized for small sec- 
ular frequencies with large electrode-ion distances. The 
latter would in addition ease the alignment of high-power 
dipole trap beams with low amount of stray light. If a 
linear RF trap is used, an improved setup shall allow to 
align the (or one of the) dipole trap beam(s) with the 
axis of the RF trap to become more insensitive to radi- 
ally defocussing DC electric fields. Alternatively, as long 
as only one or few ions have to be transferred into an 
optical trap, it might be advantageous to use an elec- 
trode geometry that provides an RF confinement in all 
three dimensions (e.g., a larger version of the stylus trap 
in Ref. 30). In such a trap DC fields would only be re- 
quired to compensate stray fields in contrast to a linear 



RF trap, where these fields also have to provide the axial 
confinement. As a result, the DC fields and thus their de- 
focussing effect could be kept weak. An additional trans- 
fer of the ion into an experimental zone could mitigate 
the problems caused by the coating of the electrodes with 
the neutral atoms from the atomic oven (cmp. Sec. I A). 
Ablation loading [31] might be an alternative. 

Displacements of the ion(s) due to light pressure can 
be avoided by using counter-propagating Doppler cool- 
ing beams, which becomes more important for further 
reduced secular frequencies. Different ion species allow 
for lower Doppler cooling temperatures, however, at the 
expense of a more complicated level scheme and thus a 
more involved optical trapping scheme. A better pre- 
cooling before transferring the ion into the dipole trap, 
e.g., by cooling it to the motional ground state [14], might 
allow for shallower optical potentials. In order to en- 
hance the lifetime of an ion in the dipole trap, cooling 
techniques used for atoms in optical traps should also 
be applicable to ions. For example, sideband or Sisy- 
phus cooling could be implemented and cavity cooling, 
as demonstrated in Ref. 32, might not even have to af- 
fect the electronic coherences. Alternatively, ions could 
be sympathetically cooled by atoms within a magneto- 
optical trap or a Bose-Einstein condensate [33] . 

We refer the reader to Ref. 1 for an extended discus- 
sions of possible applications. In the near future, cold- 
collision experiments with atoms and ions could benefit 
from the effectively micromotion-free trapping of the ion 
in an optical trap [34]. Currently, such experiments are 
implemented as hybrid setups with an ion confined in an 
RF trap and atoms in a magneto-optical trap [35] or a 
Bose-Einstein condensate [33, 36]. As the ion is subject 
to RF-driven micromotion, the system of the ion and the 
atoms experiences RF-induced heating and the collision 
energy in these experiments is limited to the order of mil- 
likelvins. Also, superimposing an optical lattice with ions 
trapped in an RF trap for quantum simulation experi- 
ments is in reach [37, 38]. On longer timescales, multiple 
ions in sparsely populated optical lattices or one (few) 
ion(s) within an optical lattice filled with neutral atoms 
might be exploited for quantum simulations. 
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